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We consider the influence of an external periodic potential on the fractional quantum Hall effect 
of two-dimensional interacting electron systems. For many electrons on a torus, we find that the 
splitting of incompressible ground state degeneracies by a weak external potential diminishes as 
cxp(— L/£) at large system size L. We present numerical results consistent with a scenario in which 
£ diverges at continuous phase transitions from fractional to integer quantum Hall states which 
occur with increasing external potential strength. 



Theoretical studies of the Quantum Hall effect (QHE) 
have identified three differentEl sources which can act sep- 
arately or in concert to create the charged excitation 
energy gap responsible for dissipationless transport in 
a two-dimensional electron system (2DES). Kinetic en- 
ergy quantization is usually primarily responsible for the 
integer QHE (IQHE), while interaction energy quanti- 
zation is responsible for the fractional QHE (FQHE). 
Integer quantum Hall gaps at fractional Landau level 
filling factors (y) can also arise from periodic external 
potentials!! Although, this instance has not yet been re- 
alized experimentally, it has played an important role 
in theoretical developments, especially in giving rise to 
the topological pictureSta of the QHE. In addition, con- 
siderations-.of the periodic external potential case have 
cautionaryu implications for the comprehensiveness of 
some theoretical pictures of the bulk QHE. 

In this paper we address the competition between 
electron-electron interactions and a periodic external po- 
tential in determining the value of the quantized Hall 
conductance of a 2DES in the strong magnetic field (B) 
limit where all electrons lie in the lowest Landau level and 
Landau level mixing can be neglected. Interactions alone 
give rise to charge gaps at a set of rational values of v, 
all of which have odd denominators, and lead to fraction- 
ally quantized Hall conductivity values <jh — {e 2 /h)v. A 
periodic external potential alone rSplits the Landau level 
into subbands in an intricate waya giving rise to the cap- 
tivating 'Hofstadter butterfly' illustrations of the mag- 
netic field dependence of the one-body spectral support. 
In the non- interacting electron limit, the charge gaps are 
the gaps in the one- bodii [Spectrum. For a given unit 
cell area, A , gaps occuirm at v = a + sA^/Aq where 
a and s are integers which depend in detail on the pe- 
riodic potential, and A^ — $o/B = 2tt£ 2 is the area 
penetrated by the magnetic flux quantum $o = hc/e. In 
this case the Hall conductivity has integer quantization, 
aii = (e 2 /h)a. Progress in nanolithography is begin- 
ning to allow the fabrication of systems in which physical 
consequences of this complicated electronic structure can 
be studied experimentally^, although the quantum Hall 
(QH) regime has not quite been reached. Addressing 
the role played by inescapable electron-electron interac- 
tions in the QH regime is a principle motivation for the 



present work.H The charge gaps of the non-interacting 
periodic potential limit clearly survive to finite inter- 
action strengths since, when all occupied subbands are 
filled, the non-interacting many-particle ground state is 
non-degenerate and separated from many-particle ex- 
cited states by a gap. The weak periodic potential limit 
is addressed below. 

The QHE of interacting electrons in a periodic exter- 
nal potential is most succinctly discussed by combining 
the toroidal geometry, in which quasiperiodic boundary 
conditions are applied to a finite area, with the. topolog- 
ical picture of the quantized Hall conductancea. In this 
picture the Hall conductance at zero temperature in e 2 jh 
units is given by the integer valued Chern index which ex- 
presses the adiabatic evolution of the phase of the ground 
state wavefunction under cyclic evolution of the bound- 
ary condition phases. The occurrence of a FQHE requires 
the many-particle ground state to be degenerate, since in 
that case the quantized Hall conductance depeads on the 
average Chern number of the degenerate statesB. Indeed, 
the necessary degeneracies are a consequenceE3 of contin- 
uous translational invariance for many-electron systems 
on a torus and occur for every state in the Hilbert space 
whether or not there is a charge gap. 

Since we restrict our attention here to magnetic field 
strengths for which charge gaps occur in both non- 
interacting modulated and interacting isotropic (i.e., zero 
periodic potential) limits, we consider only rational fill- 
ing factors v = N e /N$ — q/p where p is odd. (Here 
iV e is the number of electrons in a finite area (A) sys- 
tem and N<j, — AB/Qq is the number of states in the 
Landau level.) A gap can be created by a periodic po- 
tential alone at v = q/p only if A^/Aq = t/p for some 
integer t. A weak periodic potential will have no effect 
on the FQHE provided that the ground-state degeneracy 
splitting which it produces vanishes in the limit of an in- 
finite system size and the fractional charge gap remains 
finite. In the following paragraphs we use this criterion 
to demonstrate the stability of the FQHE against weak 
lateral potentials. 

We consider a weak periodic potential of the form 
which has been most extensively investigated for non- 
interacting electron systemsa U(r) = Uo(cos(2itx / a x ) + 
cos(27ry / a y )) . To be compatible with quasi-periodic 
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boundary conditions we choose a rectangular finite area 
system with sides L x = N x a x and L y = N y a y , where 
both N x and N y are integers, N x N y — tN e /q, and use 
a Landau gauge in which the vector potential is inde- 
pendent of the y coordinate. Our analysis is based on 
three elementary properties which follow from transla- 
tional symmetry considerations^ in the isotropic sys- 
tem limit: i) Isotropic system states can be labeled by 
a two-dimensional wavevector K in a rectangular mag- 
netic Brillouin-zone with K x 6 [0, L y /p£ 2 ) and K y £ 
[0,L x /£ 2 ); ii) Pseudomomentum is conserved so that 
(K'i'\ pp\Ki) is non zero only if K' = K + p; iii) 
Translation by (L x q/p£ 2 )y in momentum space corr 



tial QH states occurs at a critical value of Uq = U c and 
that £ diverges as U approaches U c from below. We have 
attempted to test this simplest picture for the transition 
between fractional and IQHEs, and to obtain a quanti- 
tative estimate of the modulation strength at which the 
putative transition occurs, by performing numerical ex- 
act diagonalization calculations at a series of Uq values 
including Coulombic electron-electron interactions. Here 
we report results for v = 1/3 and A^/Aq = 1/3; this 
case appears to offer the greatest promise for experimen- 
tal study since large chacg-e gaps occur in both weak and 
strong modulation limitstL3. For these commensurability 
ratios, the zero temperature Hall conductivity will change 



spond to rigid spatial translations by {2n£ 2 /L y )x so that from the fractionally quantized value a H = 1/3 e /h to 

the integer quantized value an = at the critical mod- 
ulation strength. Typical finite size results for the evo- 
lution of the low energy portion of the spectrum with 
modulation potential strength are shown in Fig. |l|. We 
expect that in the thermodynamic limit the splitting of 
the three lowest energy states should vanish and the cor- 
relation gap to the fourth lowest energy state should re- 
main finite for Uq < U c . For Uq > U c we expect the gap 
between the first and second lowest energy states, which 
is proportional to Uq for large Uq, to be finite. 



E(K,i) = E(K+L x y/p£ 2 ,i) and (K'+L x y/p£ 2 ,i'\p p \K+ 
L x y/p£ 2 ,i) — (K'i'\pp\Ki) up to a phase factor indepen- 
dent of i and i'. Here all wavevectors are understood 
to be reduced to the magnetic Brillouin-zone where nec- 
essary and pp = exp(ip ■ fj) is the density opera- 
tor. (The external potential contribution to the Hamil- 

tonian is of the form H' = {Uq/2)J2g PQ- w here G 
is a nearest neighbor reciprocal lattice vector.) When 
the FQHE occurs the p degenerate ground states (at 
K = Kq + myL x /p£ 2 with m = 0, • • • ,p — 1) are sep- 
arated from all other eigenstates by an energy-gap A 
which remains finite in the thermodynamic Umitcil. 

When H' is included, shifts in all eigenvalues occur 
starting at second order in perturbation theory. However, 
it follows from property iii) above that the shifts which 
appear at low order are identical for each member of 
the degenerate ground-state manifold. Splitting of these 
states can result only from terms in perturbation theory 
in which one member of the manifold is coupled to an- 
other via a series of intermediate states at energies above 
the gap. It follows from property ii) above that these 
terms first occur at order M y = [L x /p£ 2 ) / '(2tt / 'a y ] — 
[N x /t] for y-dependent terms in the external potential 
and at order M x = [L y /p£ 2 )/(2n/a x ] — \N y /t] for in- 
dependent terms in the external potential. (Here [i/j] 
denotes the numerator of i/j after elimination of com- 
mon divisors.) For [/o « A and fixed aspect ratio the 
splitting of the p-fold degenerate ground state manifold 
will therefore be ~ A(C/ /A) min ( M *> M v) ~ Aexp(-L/£) 
where L is the system size and £ ~ a/ln(A/£7o)- Hence, 
the FQHE will survive in an infinite system as long as 
Uq is small compared to the fractional gap A. The 
exponential dependence of splitting on system size we 
find is much weaker than the Gaussian dependence (~ 
exp(— v 2 L 2 /A£ 2 )) which occurs at the lowest ordec,pf,per- 
turbation theory for a random external potentialtlnlj and 
is in agreement with expectations from more heuristic 
argumental3 which apply equally well to random exter- 
nal potentials and are based on effective theories of the 
FQHE. 

This analysis suggests a scenario in which a continuous 
phase transition between weak and strong periodic poten- 
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FIG. 1. Low energy part of the 4-particle spectrum 
(N x = 2 = N y ) as a function of the reduced modulation 
strength V. Note the opening of the 'Hofstadter gap' Wi at 
very small modulation strength and its coexistence with the 
correlation gap W^i = Eg — Ei at a finite range of the mod- 
ulation strength which are a consequence of the finite system 
size. (Different symbols mark different pseudomomenta.) 

All the results reported here were obtained with a x — 
a y = a and up to six electrons in a system with area 



A = N X N„ 



For systems with up to five electrons the 



finite Hamiltonian matrix could be directly diagonalized, 
while for the largest systems the low lying eigenenergies 
were obtained using a block Lanczos procedure. 

In interpreting these numerical results we concentrate 
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on the energy difference W% (Wi := Ei — Eq, where Eq is 
the ground state energy) between the ground state and 
the first excited state. Explicit perturbative calculations 
for the finite size system show that at small modulation 
strengths W\ is given at leading order in perturbation 
theory by 

Wi « nN(V/A(n, N)) N A(n, N)/2 (1) 

where n — mm(N x ,N y ), N = max(N x , N y ), V — 
Uq exp(— ■nt/2p) is the external potential strength cor- 
rected by the form factor for lowest Landau level wave- 
functions, and A(n,N) = Ws(V = 0) is the finite size 
correlation gap. In the limit F » A we find that 
Wi = Ah + S(n,N) where S(n,N) is positive, and 
Ah '■= V(3 — y/S)/2 is the gap between the first and 
second subband of the Hofstadter spectrum. 
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FIG. 2. Gap between the ground state and the first ex- 
cited s tate W = W\ scaled by the dimensionless system size 
+ Ny as a function of the reduced modulation am- 
plitude v — V/A(n, N). Different curves represent different 
system sizes. The approximate crossing of these curves allows 
us to estimate the critical reduced modulation strength v c . 

We see that the finite size numerical calculations re- 
produce the expected behavior in both weak and strong 
modulation potential limits. We have performed a fi- 
nite size scaling analysisEJ to show that these results are 
consistent with a continuous quantum phase transition 
occurring between fractional {<jh = l/3e 2 //i) and in- 
teger (oh = 0) QH states at intermediate modulation 
strengths. We identify the inverse gap, W±, with the 
correlation time £r beyond which the properties of the 
system are sensitive to the periodic potential. The cor- 
relation time is finite for V$ larger than a critical value 
V c and diverges in the thermodynamic limit as Vo ap- 
proaches V c from above. W\ should obey the finite size 
scaling ansatz 



W x = {L z )-'Q{Lh^^) (2) 

where Q is the scaling function, v is the correlation length 
critical exponent, and z is the dynamical critical expo- 
nent. The small sizes of the systems for which we are able 
to numerically solve the many-electron problem limit the 
thoroughness with which we can test this ansatz. In an 
analysis which neglects aspect ratio dependence of W\ 

we use L = J N% + Ny as a measure of the system size. 

Since the fractional gap A(n, N) exhibits strong size and 
aspect ratio dependence, we perform the scaling analysis 
with the modulation strength measured in units of the 
finite-size gap defining v = Vq/ A(n, N). 
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FIG. 3. Scaling behavior of the energy gap near the criti- 
cal point. Curves for the 1 x iV-system sizes are parallel, but 
shifted with respect to each other due to next-to-leading or- 
der finite size corrections. The 2 x A^-systems show a different 
slope. Inset: Non-critical finite-size corrections in the argu- 
ment of the scaling function improve the scaling behavior. 

In Fig. H we plot L W\ as a function of v for various sys- 
tem geometries. Assumina_that z = 1 as in other quan- 
tum Hall phase transitionstj, the approximate crossing of 
these curves allows us to estimate the critical dimension- 
less modulation strength to be v c ~ 0.63. The smallness 
and limited range of our system sizes do not allow us to 
convincingly identify the correlation length critical expo- 
nent v. In Fig. 3 we show the scaling function LW\ = Q 
plotted as a function of L(v — v c )/v c , i.e. assuming v = 1. 
Deviations from scaling clearly seen in these curves be- 
come noticeably worse if v is altered by more than ~ 0.2. 
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Data points for system sizes N x x N y = 1x4 and 1x5 
fall on the same curve, consistent with the scaling ansatz. 
However, the curves for the 2x2 and 2x3, systems al- 
though parallel to each other, have a different slope than 
found for the 1 x TV systems. We can partially compen- 
sate for this large finite size effects by assuming a non- 
critical size dependence in the argument of the scaling 
function which is suggested by the small V behavior of 
Wi (see Eq. (Q)). For example, plotting Q as a function 
of n^ n 'L(v — v c )/v c , as shown in the inset of Fig. 2, 
improves the scaling behavior somewhat. 

Our perturbative expression for the correlation length 
£ suggests that the divergence of the fractional state 
correlation length £ and the vanishing of the fractional 
Hall gap (A = W3) should occur simultaneously. In our 
finite-size simulations we find that A has a minimum for 
Vq ~ V c , but at the system sizes we are able to study, no 
compelling evidence that these minimum gaps vanish in 
the thermodynamic limit emerges. 

Finally we comment on the experimental implications 
of our work. Broadly we emphasize that interactions will 
always be important at weak modulation strengths in 
any experimental investigation of a 'Hofstadter butter- 
fly' system. More particularly, we propose that continu- 
ous phase transitions between integer and FQH ground 
states will frequently occur as modulation strengths are 
weakened. Numerical exact diagonalization calculations 
for very small systems, while not able to offer compelling 
evidence, are consistent with this suggestion. For the case 
of v = 1/3 and one electron per unit cell of the periodic 
potential we estimate that the phase transition will oc- 
cur when the ratio of the modulation strength to the frac- 
tional Hall gap is ~ 0.7. T / finite-temperature scaling 
analysis of the phase transition!] implies that the Hall 
conductivity will change from an = l/3e 2 //i, to an = 
over an interval of modulation strength which vanishes 
as T~ , becoming increasingly sharp as the temperature 
is lowered. If our limited finite size data is indicative, 
vz ~ 1 rather than ~ 0.4 as found experimentally in 
field driven integer quantum Hall transitionslij. Simulta- 
neously, the longitudinal conductivity will change from 
(Ti = (corresponding to dissipationless transport when 
the fractional gap separates the ground state from the 
excited states) to a finite value at the transition point 
and back to ui = when the Hofstadter gap opens 
at modulation strengths larger than V c . These behav- 
iors would signal the continuous quantum phase transi- 
tion envisaged here. For a system with a carrier density 
n = 4 x 10 10 cm~ 2 the commensurability ratios studied 
here would be realized with a potential period a ~ 50nm 
and a field B ~ 5Tesla. We would then predict that the 
transition between integer and fractional QHEs would oc- 
cur for a modulation strength ~ O.lmeV. These periods, 
modulation strengths, and fields are not far removed from 
what can be achieved with current lithographic technol- 
ogy. The physics of the Hofstadter butterfly QH system, 
when it is finally realized, will be greatly enriched by 
electron-electron interaction effects. 
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